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Abstract 



We study the critical behavior of a general class of cubic-symmetric spin 
systems in which disorder preserves the reflection symmetry Sa — > — Sa, 
Sb Sb for b a. This includes spin models in the presence of random 
cubic-symmetric anisotropy with probability distribution vanishing outside 
the lattice axes. Using nonperturbative arguments we show the existence of 
a stable fixed point corresponding to the random-exchange Ising universality 
class. The field-theoretical renormalization- group flow is investigated in the 
framework of a fixed-dimension expansion in powers of appropriate quartic 
couplings, computing the corresponding /3-functions to five loops. This anal- 
ysis shows that the random Ising fixed point is the only stable fixed point 
that is accessible from the relevant parameter region. Therefore, if the sys- 
tem undergoes a continuous transition, it belongs to the random-exchange 
Ising universality class. The approach to the asymptotic critical behavior is 
controlled by scaling corrections with exponent A = —ar, where — —0.05 
is the specific-heat exponent of the random-exchange Ising model. 

PACS Numbers: 75.10.Nr, 75.10.Hk, 64.60.Ak 



1 



Typeset using REVTeX 



I. INTRODUCTION AND SUMMARY 



The critical behavior of magnetic systems in the presence of quenched disorder has been 
the subject of extensive theoretical and experimental study. An important class of systems 
is formed by amorphous alloys of rare earths with aspherical electron distribution and transi- 
tion metals, for instance TbFe2 and YFe2. These systems are modeled [1,2] by the Heisenberg 
model with random uniaxial single-site anisotropy, or, in short, by the random-anisotropy 
model (RAM) 



where is an M-component spin variable, Ux is a unit vector describing the local (spatially 
uncorrelated) random anisotropy, and Dq the anisotropy strength. In amorphous alloys 
the distribution of Ux is usually taken to be isotropic, since in the absence of crystalline 
order there is no preferred direction. On the other hand, in polycrystalline materials, for 
instance in the Laves-phase intermetallic (Dy-EYi_j;)Al2 compounds studied in Refs. [3,4], 
the distribution of Ux is expected to have only the lattice symmetry. 

The critical behavior, and in particular the nature of the low-temperature phase, of a 
generic system with random anisotropy depends on the probability distribution of the ran- 
dom vector Ux- In the isotropic case, i.e. when the probability distribution is uniformly 
weighted over the (M — 1) -dimensional sphere, the Imry-Ma argument [5,6] forbids the ap- 
pearance of a low-temperature phase with nonvanishing magnetization for d < A. This still 
allows the presence of a finite-temperature transition with a low-temperature phase in which 
correlation functions decay algebraically, as it happens in the two-dimensional XY model. 
Such a behavior has been predicted for the RAM with isotropic distribution in Ref. [7] and 
it has been recently supported by a 4 — e study [8,9] using the functional renormalization 
group (RG) [10]. On the other hand, standard field-theoretical perturbative approaches do 
not find any evidence for a critical behavior with long-range correlations [11-14]. While ex- 
periments have not yet found evidence of low-temperature quasi long-range order, numerical 
simulations seem to confirm the picture of Refs. [7-9], but are still contradictory as far as 
universality and behavior in the strong-anisotropy regime are concerned [15,16]. For these 
reasons, the critical behavior of the RAM can still be considered as an open problem. 

The above arguments do not apply to spin models with the discrete anisotropic distri- 
bution introduced in Ref. [11], in which the vector points only along one of the M lattice 
axes, i.e. it has the probability distribution 



where Xa is a unit vector which points in the positive a direction. This model, which we 
shall call random cubic anisotropic model (RCAM), should have a standard order-disorder 

transition: The random discrete cubic anisotropy should stabilize a low-temperature phase 
with long-range ferromagnetic order. On the basis of two-loop calculations in field-theoretical 
frameworks, it has been argued [17,18] that the transition belongs to the universality class 
of the random-exchange Ising model (REIM) for any number M of components. 
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In this paper we study the critical properties of the three-dimensional RCAM. We con- 
sider the field-theoretical approach based on the Landau- Ginzburg- Wilson kp^ Hamiltonian 

[11] 

|_ i,a ■ ijab ) 

where a,b = 1,...M and i,j = 1,...N. In the limit N ^ the Hamiltonian (3) is 
expected to describe the critical behavior of the RCAM for M-component spins. Using 
nonperturbative arguments, we show that the field theory with Hamiltonian (3) has two 
stable fixed points (FP's). One of them belongs to the REIM universality class while the 
other corresponds to the 0{N) model in the limit N 0, the so-called self- avoiding- walk 
universality class. Then, we investigate the RG flow for the model with Hamiltonian (3) 
in the framework of a fixed-dimension expansion in powers of appropriate zero-momentum 
quartic couplings. We compute the corresponding Callan-Symanzik /3-functions to five loops. 
Their analysis shows that the only accessible stable FP from the region of parameters relevant 
for the RCAM is the REIM FP. This implies that the critical behavior of the RCAM (when 
the parameters allow a continuous transition) belongs to the REIM universality class, whose 
critical exponents are [19] Ur = 0.683(3), ar = —0.049(9), rjr = 0.035(2), etc. The approach 
to the REIM scaling behavior is characterized by very slowly decaying scaling corrections 
proportional to with A = —ar ~ 0.05, which is much smaller than the scahng-correction 
exponent of the REIM, which is [20,21] A,. ^ 0.25. Our results fully confirm and put on a 
firmer ground the conclusions of Refs. [17,18] based on two- loop perturbative calculations. 

It is important to note that our results are specific of distributions that vanish everywhere 
outside the lattice axes, such as the one given in Eq. (2). Indeed, generic cubic-symmetric 
distributions P{u), and in particular the isotropic one, give rise to an additional quartic 
term that should be added to the effective Hamiltonian (3), i.e. 

Zo (t^ai<t}bi<Paj4>bj ■ (4) 
ijab 

The REIM FP is unstable with respect to this perturbation. We shall evaluate the corre- 
sponding crossover exponent, finding (f)^ — 0.79(4). Therefore, even small differences from 
the discrete distribution Pc{u) cause a crossover to a different critical behavior. Nonethe- 
less, when Pc{u) turns out to be a good effective approximation — this might be the case in 
some crystalline cubic-symmetric random-anisotropy systems — REIM critical behavior may 
be observed in a preasymptotic region. 

The general Landau-Ginzburg- Wilson Hamiltonian (3) can also be recovered by consid- 
ering systems with cubic anisotropy such that disorder preserves the symmetry s^. „ — > —s^^a, 
Sx,b Sxfi-i b ^ a. A general Hamiltonian with this property is given by 

^ - J J2 ^'-^ ■ - ^ 4,a -DoJ2J2 1x,aS%a^ (5) 

{xy) X a X a 

where — 1 and is a random vector with a probability distribution that is invariant under 
the interchange Qa <-> Qb- The exact refiection symmetry at fixed disorder — this symmetry 
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is not present in generic models of type (1) — is the key property that allows the RCAM 
and the more general class of models (5) to have a standard order-disorder transition with 
a low-temperature magnetized phase. 

The paper is organized as follows. In Sec. II we apply the replica method to the (/?^ theory 
corresponding to models (1) and (5), determining the corresponding 0^ Hamiltonians that 
are the basis of the field-theoretical approach. In Sec. Ill we discuss some general properties 
of the theory (3). In particular, we discuss the crossover behavior when randomness is 
weak, and we prove that the REIM FP is stable by evaluating its stability eigenvalues. In 
Sec. IV we investigate the RG flow by computing and analyzing the flve-loop flxed-dimension 
expansion of the /3-functions associated with the zero-momentum quartic couplings. In 
App. A we report a six-loop calculation of the RG dimensions of the bilinear operators in 
cubic-symmetric models that are used in the discussion of the stability of the FP's. App. B 
reports the proof of some identities used in the paper. 



II. EFFECTIVE *^ HAMILTONIANS 



The mapping of the RAM Hamiltonian (1) to an effective translationally-invariant 0^ 
Hamiltonian was originally discussed in Ref. [11]. In order to replace fixed- length spins 
with uncostrained variables, one performs a Hubbard- Stratonovich transformation. Then, 
for the purpose of studying the critical behavior one considers the continuum limit of the 
resulting Hamiltonian and truncates its potential to fourth order. This leads to an effective 
continuum c/?^ Hamiltonian for an M-component real field 



I 



(fx 



(6) 



where is an external spatially uncorrelated vector field with parity-symmetric distribution 
P{u) and D is proportional to Dq. We relax here the condition — 1 and require only 
that (m^) = 1, thereby fixing the normalization of D. Using the standard replica trick it 
is possible to replace the quenched average with an annealed one. The system is replaced 
by N noninteracting copies with annealed disorder. Then, by integrating over disorder, one 
obtains the following effective Hamiltonian 



repl 



/ 



(fx 



i,a ia ' ijab 



(7) 



where a, 6 = 1, . . . M, i,j = l,...N, and 

R{(f)) = -In j d^uP{u)eicp ( i:) E ^a^ft^ai^bi j ■ 



(8) 



iab 



In the limit ^ the Hamiltonian (7) is equivalent to the Hamiltonian (6) with quenched 
disorder. The expansion in powers of the field (J) can be expressed in terms of the moments 
of the distribution P{u), 



Maia2...ak = j U P{u) Ua^Ua^.-Ma^ . 



(9) 
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One obtains 
^repl = j d'^x 



tab 



iab tab ijabcd 



(10) 



Let us consider the case in which all field components become critical at Tc. This is achieved 
if the distribution P{u) is such that 



1 . 
^ M 



(11) 



This condition is satisfied if P{u) is cubic symmetric. Under this further assumption, the 
fourth moment Mabcd can be written as 



Mabcd = A {^ab^cd + ^achd + ^adh<^ + B5ahcd , 



(12) 



where A and B are parameters depending on the distribution that satisfy the Cauchy in- 
equalities A{M + 2) + B> l/M and ZA + B> It follows that 



Hrepl 



ia ia ' iab 



2M2 



[i-M^A){Y.<^ir-AD^j:. 



BD^ 



>ai<Pbi<Paj^. 



ijab 



ija 



+ 0( 



(13) 



In conclusion, for generic cubic-symmetric distributions P{u) the Hamiltonian that should 
describe the critical behavior of the corresponding RAM is 



V ia ia 



+ 



+ 4j E [(^0 + VoSij + WoSab + yoSijSab)(f)li4>lj + ZQ(t)ai(t>bi(t>aj(t>bj \ , 
ijab 



(14) 



where the term proportional to yo has been added because it is generated by RG iterations 
whenever wq ^ 0. It should be noticed that such a term arises naturally if one considers 
that, if the system is only cubic symmetric, quartic single-ion terms must be included. In 
this case it is natural to consider 



(15) 



and the corresponding (p^ Hamiltonian 



\ {d,(pf + \r(p^ - D{u . 0f + ^v,{0'f + iyo E ^'a 



(16) 
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The Hamiltonian (14) was originally introduced in Ref. [12] to describe magnetic systems 
with single-ion anisotropy and nonmagnetic impurities. 

There are two interesting particular cases. First, one may consider an 0(M) -invariant 
pure system coupled to an isotropic distribution P{u). In this case K — in Eq. (15) — 
therefore, yo — — and B — in Eq. (12), so that wq — 0. These conditions are preserved 
under renormalization by the presence of the 0{M) invariance. Note that this is not the 
case if K ^ 0, i.e. if yo ^ 0. Distributions such that B = (these distributions are not 
necessarily isotropic) give apparently Wq = 0; however, such a condition is not preserved 
under renormalization if zq ^ 0. 

A second interesting case corresponds to distributions P{u) such that A — in Eq. (12). 
It is easy to show that distributions P{u) with this property are simple generalizations of 
the distribution (2). Explicitly, they have the form 

^H = iE/K)n5K)' (17) 

where f{x) is a normalized probability distribution with unit variance. If A = 0, Eq. (13) 
implies = 0. Such a condition is stable under renormalization. Indeed, the transformation 
0ai —(f)ai for fixed a and i is a symmetry of the Hamiltonian with zq = 0, but not of the 
term proportional to zq- This symmetry is due to the fact that, for distributions of type (17), 
we can write {s ■ u)^ = I^aU^^jS^, which is symmetric under the transformations — 
at fixed u. In other words, the theory with Zq = describes models in which the refiection 
symmetry of the spins is also preserved at fixed disorder. 

In the case of discrete cubic-symmetric distributions of type (17), we have 

^0^^, Wo^-12BD'. (18) 
Apparently, these conditions imply 

1*0 > 0, Wo< 0, Muo + Wo < 0, (19) 

where the last condition follows from the bound B > 1/M. The equality Muq + wq — 
is obtained by using distribution (2). Relations (18) and (19) should be considered as 
indicative, since the mapping between the Hamiltonian (16) and the general Hamiltonian 
(14) gives also rise to higher-order terms. 

It is also interesting to consider the effective continuum Hamiltonian corresponding to 
Eq. (5). In this case we obtain 



(20) 



If P{q) is invariant under the interchanges qa ^ Qb, we can write for the first moments 
Ma = a and Mab = b + c6ab- A simple calculation gives again the general Hamiltonian (14) 
with Zq — and 

u^^—{a^-h), ^0 = -^- (21) 
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Since b + c> a^, we obtain 

uo + Wo < 0. (22) 

Equality is obtained for P{q) = Yia^i^a — 1) (in this case however Wq = 0). Finally, note 
that if c = then we have wq = 0. Such a condition is stable under renormalization, and 
thus this class of models is expected to have a different critical behavior. It corresponds to 
the one of the randomly dilute cubic models discussed in Ref . [22] . Distributions with this 
property are however quite peculiar. They have the general form 

M 

^(?)=/(?i)n (23) 

a=2 

The stability region of the quartic potential in the </7^ Hamiltonian (3) is given by the 
conditions 

Nuo + vo + Nwo + yo>0, (24) 

NMuo + Mvo + Nwo + yo>0, (25) 

uo + vo + wo + yo> 0, (26) 

Muo + Mvo + wo + yo>0. (27) 

However, as discussed in Ref. [12], in the zcro-rcplica limit N ^ 0, the only relevant stability 
conditions are those obtained by using replica- symmetric configurations. Therefore, for the 
RCAM one should only consider Eqs. (24) and (25) with = 0, i.e. 

Vo + yo>0 if t'o > 0, 
Mvo + yo>0 if < 0. (28) 

Equivalently, the relevant stability conditions can be obtained by considering the Hamilto- 
nians (16) and/or (20). 



III. GENERAL RENORMALIZATION- GROUP PROPERTIES 

A. Fixed points of the theory 

The properties of the RG flow arc essentially determined by its FP's. Most of them 
can be identified by considering the theories obtained when some of the quartic parameters 
vanish. For example, we can easily recognize: 

(a) the 0(M x N) theory for Vq = Wq = yo = 0; 

(b) N decoupled 0(M) theories for uq — Wq — yo — 0; 

(c) M decoupled 0{N) theories for uq — Vq — yo — 0; 

(d) M X decoupled Ising theories for uq = vq = wq = 0; 

(e) the MN model (see, e.g., Refs. [23,24]) for wq = t/o = 0; 

(f) the NM model for vq = yo = 0; 

(h) N decoupled M-component cubic models for uq — Wq — 0; 
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TABLE I. Fixed points of the Hamiltonian (3) near four dimensions. We report the lead- 
ing nontrivial contribution of the expansion in powers of e, taken from Refs. [11,17]. Here, 
Kd = (47r)'^r(d/2)/2, a± = (M - 4 ± + 48)/8, /3± = -(M + 12 ± VM2 + 48)/6, 

A±± = 6a± + 3/3± + M + 6. The general expressions for the stability eigenvalues of FP's XI-XIV 
are rather cumbersome. We only report their numerical value for M = 3. 



y* /K^ Stability eigenvalues 



I 


Gaussian 





II 


0(M) 


M+8'= 


III 


0(0) 





IV 


0(0) 





V 


Ising 





VI 




3 




2(M-1) 


VII 







VIII 


Cubic 


2 

m'- 


IX 


M ^ 2 




X 







XI 






XII 






XIII 






XIV 






XV 


REIM 





XVI 


REIM 





XVII 


M = 2 








h 





3(M-4) 
8(M-1) 

h 



M-4 
4(M-2) ^ 

M 

i^" 
3tx+ 

i^" 
3a_ 

3a_ 















-h 

3{M+4) 

4A+ + 
3(M+4) 

4A+_ 
3(M+4) 

4A_+ 
3(M+4) , 

4A ^ 











h 




2(M-4) 
3M *= 
M-4 
3{M-2) ^ 

¥ 

A- + 

4^: 

-1-4 , /U /; 
=■=3 V 53 ve 



UJu - 
UJu ■ 

Ulu ■ 

■ 



_ A-M _ 4+M 

e,u}u — —jY^:ge,uiw — - jgqige,'^!/ 
■- e,u)v = uiw = uiy = e/2 
: Wv = — €/2,ajui = e,Ljy = e/2 
: oJi, = LJu, = — e/3, ojy = e 
4-M 

€, 0)2 — tDj, — jpj^rrje, u)u) 
o)^ = e, = 

. ... 4+A<f 
0)2 — Tn~^''^w — 3M^^' 4 

' (4+M) 



M-H8'' 



4+M 
4(M-1) ' 



4-M 
3M 
_ 4-M 
^ ^'^^2 - 6(M-2)"=' 



_ -(4-M) 
(M-2) ' ^4 - 6(M-2) ' 



u>y = -e 

e,u>m = 

e, ijJv =ijJ3 = e/3, uji = —e/3 
e,W2 = 1.33e,u)3 = 0)4 = -1.43e (for M = 3) 
e,a)2 = -^3 = 0.371e,a)4 = -0.344e (for M = 3) 
e,a;2 = -0)3 = 0.435, W4 = -0.403e (for M = 3) 
e,W2 = 0)3 = -3.32e,a;4 = -3.08e (for M = 3) 



: w„ = = i-y/H'ye, 0^2 = 2e 
= a;„ = -wi = iFy^ye, u)2 = 2e 



: 013 = 



53 ' 

/"2T ^ 



2e 



(i) M decoupled A^-component cubic models for uq = Vq — 0; 
(j) M X iV-component cubic model for = = 0; 

(k) the randomly dilute M-component cubic model (see Ref. [22]) for wq — Q and N — 
(1) the tetragonal model [23,25] for M = 2 and two = 0. 

The FP's of these theories arc also FP's of the enlarged model (3). Of course, there 
may also be FP's that arc not related to the above particular cases. Their presence can 
be investigated by low-order e expansion calculations. First-order e-expansion calculations 
[11] show the presence of 14 FP's for M 7^ 2 and of 13 FP's for M = 2. As in the 
REIM case, at two-loop order other 0{^/e) FP's appear [12]: 4 FP's for M 7^ 2 and 6 
FP's for M = 2 [17]. The two-loop e-expansion results are summarized in Refs. [17,18]. 
In Table I we report the leading e-expansion terms for the location of the FP's (in the 
minimal-subtraction renormalization scheme) and the corresponding stability eigenvalues. 
The only stable FP's are the 0(0) and the REIM FP's, which are already present in models 
(a) and (i), respectively. These results have also been supported by two- loop fixed-dimension 
calculations [17]. In order to understand the relevance of the various FP's for the RCAM, 
we need to check which one is accessible from the region of the quartic parameters relevant 
for the three-dimensional RCAM. This issue will be investigated in Sec. IV computing and 
analyzing five-loop series in the framework of the fixed-dimension expansion. 
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B. Crossover behavior close to the pure spin model 



The 0(M) -symmetric FP located in the v-a.xis describes the critical properties of the pure 
spin system in the absence of cubic anisotropy. It is interesting to compute the crossover 
exponent in the presence of random anisotropy. Setting tp = {T — Tp)/Tp, where Tp = 
Tc{Dq = 0) is the critical temperature in the absence of anisotropy, in the limit tp ^ and 
Dq ^ the singular part of the free energy can be written as 

J^^I^.I^-^/Po^tl"^"), (29) 

where Ut ^ tp + aiD^ + a2tp + . . . is the scaling field associated with temperature, a is the 
specific-heat exponent in the 0(M) theory, is the crossover exponent, and f{x) is a scaling 
function. As a consequence of Eq. (29), for sufficiently small Do the critical-temperature 
shift is given by 

ATe(Do) = T,{Do) - Te(0) ^ aDo'/'^^ +bDl + cD^ + ... (30) 

The crossover exponent do is related to the largest positive RG dimension of the perturba- 
tions at the 0(M) FP that are present in the Hamiltonian (3), i.e. of the terms proportional 
to Mo, Wq, and yg. For uq = Wq = yo = 0, the Hamiltonian (3) describes N decoupled 
0(M)-symmetric systems. The RG dimension of the terms proportional to uq and wq can 
be determined by writing [11] 

^(tiO + WQ5ab)(t>li(t>lj = M{Muo + Wq) ^ £i8j + Wo X! TaaiXaaj, (31) 
abij ij ija 

where 

^« = T7 X] ^ai^ 

a 

Tabi = 4>ai4>bi — ^ab^i- (32) 

The bilinears £i and Tau are respectively the energies and the quadratic spin-2 operators of 
the N decoupled models. \i Ve = ^/^ and y^ are the corresponding RG dimensions, the 
two terms given above have RG dimensions yu = ^yE — 3 = a/u and y.^, = 2y'j' — 3. The 
perturbation proportional to yo does not couple the different replicas and therefore its RG 
dimension yy is simply the RG dimension of the cubic perturbation, which is related to the 
RG dimension of the spin-4 perturbation of the 0{M) FP [26,27]. Therefore, the 0(M) 
FP is perturbed by three terms of RG dimensions y^, yw, and yy, which can be determined 
using known results for the RG dimensions of generic operators in an 0(M) theory, see, 
e.g., Refs. [28,25] for reviews of results. Since a is negative for M > 2, we have |/„ < and 
therefore the corresponding term is always irrelevant. The exponent yx has been obtained 
by using field-theoretical [26] and Monte Carlo methods [29]: field-theoretical analyses give 
yT = 1.766(6) for M = 2 and yr = 1.790(3) for M = 3, while Monte Carlo simulations 
give yT = 1.756(2) for M = 2 and yx = 1.787(2) for M = 3. Correspondingly, we find 
yyj = 0.532(12) and = 0.511(6) for M = 2, and y^j = 0.580(6) and 2/^ = 0.573(3) for 
M — 3. Therefore, the perturbation proportional to wq is always relevant. Finally, using the 
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results of Rcfs. [26,27] for the spin-4 perturbations at the 0(M) FP, we have Uy = —0.103(8) 
for M = 2 and yy = 0.013(6) for M = 3. This impHcs that the y-term is irrelevant for M = 2, 
but relevant for M = 3. In conclusion, for both M = 2 and M = 3, the most relevant quartic 
perturbation is given by the w-term, which determines the crossover from the pure critical 
behavior in the limit of small anisotropy strength. Therefore, (po — Vw^ — 0.357(3) for 
M = 2 and (po = Vwi^ = 0.412(3) for M — 3. In the crossover limit in which Dl\ut\~'^'^ is 
held fixed, the operators with RG dimensions yu and yy give rise to scaling corrections. In 
particular, there are corrections proportional to f^^ , with Ay = yyU — 4>d, Ay = 0.426(6) 
for M = 2 and A^^ = 0.403(5) for M — 3, which are more important than the usual 0{M)- 
invariant corrections, which vanish as t^, with [25] A fa 0.54 for M = 2 and A !v 0.56 for 
M = 3. 

It is worth mentioning that the scaling behavior (29) with the same crossover exponent 
(pD also holds for a RAM with generic distribution P{u), and in particular for the isotropic 
case. Indeed, the additional term proportional to Zq appearing in the Hamiltonian (14) has 
the same RG dimension of the w-term at the 0(M) FP. This can be inferred by rewriting 

Yl (Pai<f>bi<f>aj4>bj = Yl '^abi%bj +^12 ^^^J' (3^) 
abij abij ij 

where Tau and Si are defined in Eq. (32). The first term is the most relevant one and 
therefore, we obtain y^ = 2yT — 3 and also yz = Vw 

Let us note that the relatively small value of (f)D makes the measurement of (f)D from the 
critical-temperature shift for small random anisotropy rather difficult. Indeed, in Eq. (30) 
the term D^'^" is suppressed with respect to the first two analytic terms proportional to 
Dl and D^, since 2/0/) ^ 4.9 (resp. 2/0^ 5.6) for M = 3 (resp. M = 2). This explains 
the results of Ref. [4] that measured in crystalline Laves-phase (Dya.Yi_2,)Al2 for different 
values of x. Since [30] Z^o as x ^ 1, they were able to measure ATc{Dq) for — > 0. 
The experimental results were fitted assuming 

ATe(Do) ~ dI^^: obtaining i/j = 0.80(8). 
This result is in substantial agreement with the theoretical prediction ATc{Dq) ~ Dq, but 
does not provide information on the crossover exponent 0d- 

For M > 3, pure systems with Hamiltonian (15) do not have a critical behavior in the 
0{M) universality class, see, e.g., Refs. [23,27,22]. If the system has [111] as easy direction, 
its critical behavior belongs to a different universality class with reduced cubic symmetry, 
while systems with [100] easy axis are expected to show a first-order transition. In the 
latter type of systems, randomness may soften the first-order transition. This issue shall be 
discussed in Sec. IV C. On the other hand, we now show that in cubic systems with [111] 
easy axis randomness is a relevant perturbation and therefore, for small randomness, these 
systems show a crossover behavior with positive exponent 0d, cf. Eq. (29). For uq = Wq — 
the Hamiltonian (3) reduces to the one for N decoupled systems with cubic symmetry. The 
RG dimensions of the terms proportional to uq and Wq at the cubic FP provide the crossover 
exponent (f)^- In order to determine them, we use again Eq. (31). The RG dimension of 
£i is yE = l/i^, where u is the correlation-length exponent, while that of Uai = laai, Uu, 
is computed in App. A by resumming its six-loop perturbative expansion. Thus, the RG 
dimensions of the two terms appearing in the right-hand side of Eq. (31), we denote them by 
yu and y^, are given by y^ — 2yE — 3 — a / 1/ and y^; — 2yu — 3, respectively. Since a < at 
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the cubic FP for any M > 3, the first term is irrelevant. On the other hand, the estimates 
of uu reported in App. A show that > for any M > 3. For example — 0.549(14) for 
M = 3, and therefore (po = 0.387(14). 

Note that in a generic cubic-symmetric RAM, one should also consider perturbations 
proportional to zq. We use again Eq. (33). However, in the presence of cubic symmetry Tabi, 
cf. Eq (32), is not an irreducible tensor. One must consider separately Uai = Taai and Tabi 
with a ^ b, that have different RG dimensions yu and yT- Therefore, the term proportional 
to Zq is the sum of three terms of RG dimensions 2yE — 3 = a/i', 2yu — 3, and 2yT — 3. The 
last one is the largest, so that yz — 2yT — 3. Using the results of App. A for M = 3, we 
find yz — 0.600(4). The exponent is larger than yyj. Therefore, a generic cubic-symmetric 
RAM shows a different crossover behavior with crossover exponent (f)D — yz^ — 0.427(3). 

C. Stable fixed points 

The critical behavior in the presence of random anisotropy should be described by the 
stable FP of the theory (3) which is accessible from the RGAM. The two-loop e-expansion 
calculations of Ref. [17] summarized in Sec. Ill A find two stable FP's. One of them is 
located in the w-axis, and it is associated with the 0(0) or self- avoiding- walk universality 
class. This FP is also stable in three dimensions. Indeed, the terms proportional to Wq, Wq, 
and ai'G interactions transforming as the spin-4 representation of the 0(M x N) group. 
Therefore, they have the same RG dimension which is given by yv,w,y = —0.37(5), obtained 
in Ref. [26] from the analysis of six-loop fixed-dimension and five-loop e series. It was argued 
in Ref. [17], on the basis of two-loop calculations, that the 0(0) FP is not accessible from the 
parameter region relevant for the RGAM. This will be confirmed by the five-loop analysis 
of the RG flow presented in Sec. IV. 

For Mo = Wo = the Hamiltonian (3) corresponds to a cubic-symmetric model and, for 
A?" — > 0, it is the sum of M independent models that are the field-theoretical analog of the 
REIM. We will now show that the REIM FP is stable in the theory (3). It is sufficient to 
show that the terms proportional to Uq and Vq are irrelevant. For this purpose, we rewrite 

^(«o + vo5,,)<t^l^l = N{Nuo + t;o) E^a + E^a., (34) 

ijab 0, ai 

where Sa = jjHi <i>i,i and Uai = 4>li~^a- For N ^ 0, £i and Uai have the same RG dimension 
(see App. A for the proof), ys — yu — ^/^r, where Ur is the correlation-length critical 
exponent of the REIM universality class. Therefore, the RG dimension of the perturbation 
is given by yuv = 2yE — 3 = ar/i^r, where a,, is the REIM speciflc-heat exponent. Since ar 
is negative, see the estimates reported in Refs. [25,31,19], the REIM FP is stable. Using the 
recent Monte Garlo results reported in Ref. [19], we finally arrive at the estimate y^v ~ —0.07. 
As we shall see in Sec. IV, the REIM FP turns out to be accessible to the RGAM, and no 
other stable FP exists in the region relevant for the RGAM. Therefore, the REIM universality 
class describes the critical properties of the RGAM in the case it undergoes a continuous 
transition. Estimates of several universal quantities for the REIM universality class can 
be found in Refs. [25,31,19,32]. Note, however, that the critical exponent controlling the 
leading scaling corrections differs from the one for the REIM, which is [20,21] ~ 0.25. In 
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the RCAM the leading scahng correction is due to the Hamiltonian terms proportional to 
Uq and Vq. They cause slowly decaying corrections of order with 

A = -a^ = 0.049(9). (35) 

If P{q) = HaPaiQa), i-^. the probability distributions of the variables qa are independent, 
the stability of the REIM FP can also be proved by starting directly from Eq. (20). Indeed, 
such a Hamiltonian corresponds to M random-exchange ip'^ models coupled by the term 
proportional to vq. Such a term has the form Y^ab^a^b, where £a — jf(fia is the energy 
of the REIM. Therefore, this perturbation has RG dimension 2/z/f. — 3 = ar/^r-, which is 
negative. Thus, the coupling among the models is irrelevant, and thus it does not change 
the universality class of the system. If P{q) does not factorize, the M (p^ models are also 
coupled by disorder. The above-reported analysis shows that also this coupling is irrelevant, 
its RG dimension being ar/^r < 0- 

As discussed in Sec. II, in the case of a generic random cubic-symmetric distribution 
P{u), the Hamiltonian (14) also contains the term proportional to z^. It is important to 
note that the REIM FP is unstable with respect to this perturbation, since its RG dimension 
Hz is positive at the REIM FP. The dimension of this perturbation, y^, can be computed by 
rewriting the term proportional to zq as 

XI 4>ai4>bi4>aj4>bj = X X '^o.ijTbij + ^aiUbi + SaSb, (36) 

abij ab i^j ab i ab 

where Taij = 4>ai4>aj with i 7^ j. Therefore, this perturbation is the sum of three terms that 
have RG dimensions 2yT — 3, 2yu — 3, and 27/^ — 3. Using the results reported in App. A, 
one finds that the most relevant term is the first one, so that 

= 2yT - 3. (37) 

Using the estimate yr — 2.08(3) reported in App. A, one obtains 

= 1.16(6), 0, = |/,i/ = 0.79(4), (38) 

where 0^ is the corresponding crossover exponent. 

D. Critical behavior for infinitely strong random anisotropy 

In this Section, we wish to investigate the general model (5) in the limit of infinite 
disorder, showing that, under some mild hypotheses for the probability distribution P{q), 
one has REIM critical behavior for M = 2 and M — 3 and no transition for M > 4. This 
analysis further confirms the results of Sec. HI C. 

We first consider the case Dq — * +00. We suppose that the distribution P{q) is such 
that there is only one direction k such that q^ = max^, q^ (or at least that this condition 
is verified with probability one). This is the case if the distribution is continuous and is 
also true for the distribution P{q) derived from (2) (note that qa = ul). Because of the 
assumption on P{q), for Dq — > -|-oo the spin s is constrained to lie along the k direction. 
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i.e. Sk — ±1, Sq = for a ^ k. Thus, in this hmit we can rewrite the Hamiltonian in the 
following way. At each site we define M Ising variables a and M disorder variables p^^a- 
The Ising variables assume values ±1, while the disorder variables assume values and 1 
with probabilities induced by the distribution of q: 

Px,a = 1 if qx,a > qx,b for every 6 7^ a , 
Px,a — otherwise. 

Then, the average value of a quantity 0{sx,a) is given by 

JO) = [{Oip,,aa.,a))al, (39) 

where [-j^ indicates the average over the disorder variables px,a and indicates the sample 
average with Hamiltonian 

- ^x,a(^y,aPx,aPy,a ■ (40) 

a {xy) 

If O depends only on a single component of the spins, say it depends only on s^^i, we can 
integrate out a^^a a-nd Px,a for a > 2. Thus, the Hamiltonian becomes a REIM Hamiltonian 
with disorder p^-j. Now, we use the symmetry of P{q) to conclude that the probability 
that Pa is 1 must be independent of a. Since "^^aPa is always equal to 1, we obtain that 
that px,i = 1 with probability 1/M and px,i = with probability 1 — 1/M. Therefore, we 
obtain that correlation functions of si are exactly equal to the correlation functions of the 
site-diluted Ising model with vacancy density 1 — 1/M. Note that this result is not true for 
correlation functions that involve different components of the spins. Indeed, the Hamiltonian 
(40) corresponds to M REIM, but they are coupled by the disorder variables. Thus, these 
correlation functions are not simply obtained by multiplying REIM correlation functions. 
These considerations allow us to predict the behavior of the model (5) for Dq — > +00. Since 
the REIM has a continuous transition for spin density p > Pc, Pc = 0.3116081(13) on a cubic 
lattice [33], we predict that the model has a continuous transition for M — 2 and M — 3 
and no transition at all for M > 4. 

Let us now consider the opposite case Dq — —00. If the distribution P{q) is such that 
there is only one direction k such that qk — min^ qi,, the previous argument applies with no 
changes. Note that distribution (2) does not satisfy this condition for M > 3. Indeed, in 
the limit Dq — ^ —00 the spins are constrained to be orthogonal to q, and therefore cannot 
be considered as Ising variables. In this particular case, the behavior at the transition, if it 
exists, is not predicted by this argument. 

IV. RENORMALIZATION-GROUP FLOW IN THE QUARTIC-COUPLING 

SPACE 

A. The fixed-dimension five-loop expansion 

In this section we study the RG flow of the theory (3), determining the stable FP's and 
their attraction domain. For this purpose, we determine the five-loop perturbative expansion 
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of the j3 functions in terms of appropriately defined zero-momentum quartic couplings at 
fixed dimension. In the present case we define v, w, and y by writing 



^aibjckdli^) — ^^4> RMN^aibjckdl + ^ RMBaibjckdl + ^ RNCaibjckdl + U ^aibjckdl) , (41) 




where Rk — 9/{S + K), A, B, C, and D are appropriate tensors defined so that at tree 
level, uq — muRMN, Vq = itivRm, wq = mwRN, and j/o — "f^y- The mass m and the 
renormalization constant are defined by 



Here T^^) and T^^) are the four- and two-point one-particle irreducible correlation functions. 

We computed the /3-functions to five loops, which required the calculation of 161 Feyn- 
man diagrams. We employed a symbolic manipulation program, which generated the di- 
agrams and computed the symmetry and group factors of each of them. We used the 
numerical results compiled in Ref. [34] for the integrals associated with each diagram. We 
do not report the series for generic values of and M, but only for the the physically inter- 
esting cases = and M = 3, 2. The series for generic N and M are available on request. 
The coefficients of the five- loop series of the /5-functions and Py are reported in 

Tables II- V for M = 3 and Tables VI-IX for M = 2. We performed the following checks: 

(i) The series are symmetric under the transformation v ^ w and M ^ N; 

(ii) For w — y = and for w = y = (with N <-> M), we recover the series for the MN 
model reported to six loops in Ref. [24]; 

(iii) For v — w — we obtain the series for the (M x A'")-component cubic model, reported 
to six loops in Ref. [27] ; 

(iv) We obtain the series for the randomly dilute cubic model ior w — and N — and 
those of the tetragonal model for ti; = and M — 2. These series are reported to six 
loops in Refs. [22,25]. 

(v) For N — the series satisfy the identities 



where Preim,u{u, y) and /9reim,2/(w, y) are the /3-functions of the REIM model obtained 
by setting v — w — Q and M — 1. The corresponding six- loop series are reported in 
Ref. [27]. These identities are proved in App. B. 

(vi) At two loops the series agree with the expansions reported in Ref. [17]. 

Perturbative series are divergent and thus a careful analysis is needed in order to obtain 
quantitative predictions. In the case of systems without randomness, they are conjectured 
to be Borel summable and this allows one to use the Pade-Borel method or methods based 



(42) 



(3u{u, 0, w, y) + (3w{u, 0, w, y) = /3reim,«(h + w,y), 

Py{u, 0, W, y) = pRElM,y{u + w,y), 



(43) 
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on a conformal mapping [35]. In random systems, the pcrturbative approach faces addi- 
tional difficuhies: the pcrturbative series are expected not to be Borel summablc [36,37]. 
Nonetheless, in the REIM case quite reasonable estimates of the critical exponents have been 
obtained by using the fixed-dimension expansion in d = 3 (see, e.g., Refs. [25,31]). Similarly, 
the usual resummation methods applied to the RCAM expansions give quite stable results, 
at least when the quartic couplings are not too large, giving us confidence on the correctness 
of the conclusions. 



B. The RG trajectories 

The knowledge of the (3 functions allows us to study the RG flow in the space of the 
quartic renormalized couplings u, v, w, and y. For this purpose we follow closely Ref. [21]. 
The RG trajectories are lines starting from the Gaussian FP (located atu — v = w = y = 0) 
along which the quartic Hamiltonian parameters Uq, Vq, Wq, and yQ are kept flxed. The RG 
curves in the coupling space depend on three independent ratios of the quartic couplings. 
The RG trajectories can be determined by solving the differential equations 

-A^^A., (44) 
where gi — u, v, w, y, and A e (0, oo), with the initial conditions 

= Si, (45) 

A=0 

where si = s„ = uo/\vq\, S3 = = w;o/|^'o|, S4 = Sy = yo/\vo\, and S2 = +1 if vq > 0, 
S2 = —1 if Vq < 0. The functions gi{X, Si) provide the RG trajectories in the renormalized- 
coupling space. The attraction domain of a FP g* is given by the values of uq, vq, wq, and yo 
corresponding to trajectories ending at g*, i.e. trajectories for which gi{X = 00, Si) = g*. We 
recall that the 0(M) FP is located in the v-axis at v* ^ 1.40, 1.39 for M = 2, 3 respectively 
[38,39]; the 0(0) FP lies in the w-axis at u* ^ 1.39 (Refs. [38,39]); the REIM FP lies in 
the w-y plane at w* ~ —0.7 and y* ~ 2.3 (we report here the field-theoretical estimates of 
Ref. [24]; Monte Carlo estimates are given in Ref. [19]). 



C. Results 

In this Section we report our analyses of the five-loop pcrturbative scries. We have 
resummed the /^-functions by using the Pade-Borel method. The major numerical problem 
we faced was the fact that most of the approximants were defective in some region of the 
couphng space, forbidding a complete study of the RG fiow. This is not unexpected since 
the pcrturbative series are not Borel summable. Approximant [3/1] for P^, [4/1] for and 
[3/2] for j3w with 6 = 1 [40] were not defective in all the region of the RG flow we considered 
(in some cases the Pade approximant to the Borel transform had a pole on the positive real 
axis but far from the origin, in a region that gives a negligible contribution to the resummed 
function). On the other hand, all approximants for (3y were defective somewhere in the 
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FIG. 1. Projections of the RG flow for the three-component case, M = 3, in the y-w, y-v 
and y-u planes, as a function of Su = uq/vq, for uq > 0, vq > 0, and wq < 0. Here yo = and 
Suo + wo = 0. The REIM FP corresponds to u* = v* = 0, w* -0.7, and y* 2.3 (Ref. [24]). 

region we wished to investigate. For jSy we used approximant [3/2] with 6 = 1, that had the 
least extended defective region. All results we present here were obtained by using these 
approximant s. It must be stressed that other choices gave results that were similar in the 
regions in which they were well-defined. 

First, we checked the general results reported in Sec. III. We considered the 0(M), cubic, 
0(0), and REIM FP's and for each of them we determined the stability eigenvalues. The 
results are in full agreement with the conclusions of Sec. Ill, confirming that the 0(0) and 
the REIM FP's are stable. Then, we looked for additional FP's beside those identified by 
the e-expansion analysis of Sec. Ill A. For this purpose we considered the RG flow starting 
from arbitrary values of u,v,w,y. We only observed runaways trajectories or a flow towards 
either the REIM or the 0(0) FP's, confirming that the REIM and the 0(0) are the only 
stable FP's. In particular, trajectories corresponding to Hamiltonian parameters wq < 0, 
Mo > 0, uq + Wq < 0, and that satisfy the stability bound (28) never flow towards the 0(0) 
FP, which is therefore not accessible from this region. They either flow towards the REIM 
FP or apparently run away towards inflnity. 

For the purpose of illustration, we flrst consider the case yo = 0, Wq/uq = —M, and 
Vq > 0, which apparently corresponds to the model (6) with distribution (2), cf. Eq. (18) 
(note that B = 1/M in this case). In Fig. 1 we show the RG trajectories for M = 3 for 
several values of Su > 0. The approximant of Py is defective for 0.05 ^ Su ^ 0.3 and y close to 
1. This explains the sudden change of direction of the trajectory with Su = 0.3 in Fig. 1 when 
y is close to 1. For M = 3 (resp. M = 2) the RG trajectories appear to approach the REIM 
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FP for < s„ ^ 0.9 (resp. < ^ 1.4). For larger vahics of Su the flow runs close to regions 
in which some approximant is defective. Apparently, trajectories flow towards infinity, but 
this could be an artifact of the resummation. In any case, if true, this would imply that the 
corresponding systems do not undergo a continuous transition. As a consequence, since Su is 
directly related to the anisotropy strength D, the continuous transition would be expected 
to disappear for sufficiently large values of D. These conclusions do not immediately apply 
to fixed-length spin systems, i.e. to the Hamiltonian (1) since in this case [41] vq = +oo. 
Thus, it is not clear which is the correct value of even for uq = oo. The critical behavior 
of this system for strong disorder has been discussed in Sec. Ill D. 

The qualitative picture does not change if we do not require Wq/uq — —M and Sy — 0. 
For instance, we can consider the case wq/uq — —M, vq > 0, and arbitrary Sy. We are 
able to resum reliably the perturbative series for Sy > —0.7 and there we observe that some 
trajectories flow towards the REIM FP, while others run away to infinity. The attraction 
domain of the REIM FP enlarges with increasing Sy-. it is approximarely bounded by < 
0.9 + O.Gsy for M = 3 (s„ < 1.4 + lAsy for M = 2) in the region -0.6 < Sy < 0.3 
(—0.7 ^ Sy ^ 1). For larger value of Sy, the attraction domain becomes even larger and 
extends beyond the lines reported above. For Sy ^ —0.6 some approximants become defective 
and we cannot determine reliably the RG fiow. As in the case Sy = 0, there is some evidence 
that trajectories flow towards inflnity for s^^ ^ — 1, while for — 1 < < — 0.6 they may still 
flow to the REIM FP. 

Then, we have investigated the behavior for vq < 0, although this region does not appear 
to be of physical interest. As in the pure case, all trajectories apparently run away to infinity. 

Finally, let us discuss whether 0(0) critical behavior can be observed by appropriately 
tuning the model parameters. We have investigated this question in detail. We find that 
the 0(0) FP can be reached only if uq + Wq > 0, irrespective of the other parameters as 
long as Mo > and Wq < 0. This result can be proved straightforwardly in the limiting case 
Vq = 0. Indeed, since Z^^, = for f = 0, if we start with vq = the flow will be conflned 
in the hyperplane v = 0. But for f = we can use identities (43) that show that the flow 
for the couplings u + w and y is identical to the flow observed in the random-exchange (p'^ 
theory. Therefore, for uq + Wq — we observe pure Ising behavior, while for uq + Wq < 
(resp. Mo -l- Wq > 0) RG trajectories flow towards the REIM (resp. 0(0)) FP. Therefore, 
if Wq/uq < —1, as implied by Eq. (22), only REIM critical behavior can be observed. For 
Vq > 0, similar conclusions are obtained numerically: The attraction domain of the 0(0) FP 
is included in the region Wq/uq > — c, where the constant c is positive and smaller than 1, 
depends on s^, and tends to 1 as — > — oo, i.e. Vq — > 0. 

In conclusion, our analysis gives a full picture of the critical behavior for cubic magnets 
that have vq > — we have vq = +oo for flxed-length spins [41]. For M = 2 the pure system 
has a critical XY transition for Sy > —2/3, an Ising transition for Sy = —2/3, and a flrst- 
order transition for —1 < Sy < —2/3 [values of the parameters such that Sy < —1 are not 
allowed since they do not satisfy the stability bound (28)]. Note that first-order transitions 
cannot be observed in the pure model (15) with fixed-length spins. Indeed, for the strongest 
possible negative anisotropy, K = — oo, the Hamiltonian can be written as two decoupled 
Ising models [42], and thus the system with K = — oo exactly corresponds to Sy = —2/3. 
As a consequence, finite values of K have Sy > —2/3, and therefore the model is expected 
to have always an XY transition. Randomness changes the critical behavior. For small 
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randomness and small anisotropy, we always predict REIM critical behavior, while for large 
disorder (unless we consider fixed-length spins, cf. Sec. HID) we do not expect a continuous 
transition. Note that the behavior of systems with —l<Sy< —2/3 remains unclear since 
in this region we are not able to resum reliably the perturbative expansions. In particular, 
we cannot clarify if softening occurs. A Monte Carlo simulation [43] found that model (15) 
with K = — oo has a continuous transition for small disorder, in agreement with our results. 

For M = 3 we expect a continuous transition for Sy > and a first-order one for Sy < 0. 
If we add randomness to systems with Sy > 0, the continuous transition survives but now 
belongs to the REIM universality class; for large disorder the transition may disappear. For 
Sy < one may observe softening, i.e. the first-order transition may be changed into a 
continuous one by small disorder. Note that softening always occurs for infinite disorder, 
Dq = +00, in the model (5), independently of the sign of yo, under mild assumptions on the 
distribution P{q), see Sec. HID. 
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APPENDIX A: RENORMALIZATION-GROUP DIMENSIONS OF BILINEAR 
OPERATORS IN THE CUBIC-SYMMETRIC THEORY 



In this appendix we compute the RG dimensions of bilinear operators in the cubic- 
symmetric theory 



He 



(fx 



1 

4! 



TV 



(Al) 



where </? is an A'"-component field. The bilinear operators can be written in terms of tensors 
belonging to different irreducible representations of the cubic group: 
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TV 
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(A2) 

(A3) 
(A4) 



The RG dimension of the energy operator is = 1/z/, where u is the correlation-length 
exponent. The RG dimensions of the operators Ui and Tij, respectively yu and yr, in the 
cubic-symmetric theory (Al) will be computed below. Note that in O(A^) -symmetric theories 
the tensors Ui and Tij belong to the same irreducible representation and therefore yr — yu- 
In cubic systems this is no longer the case. 

In order to compute yu and yx, we consider the perturbative approach in terms of the 
zero-momentum quartic couplings u and v at fixed dimension. We refer the reader to Ref. [27] 
for notations and definitions; there one can also find the six-loop perturbative expansion of 
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the /?-functions and of the RG functions associated with the standard exponents. In order to 
compute the RG dimensions of Ui and T^, we consider the related RG functions Zjj and Z^, 
defined in terms of the zero-momentum one-particle irreducible two-point functions Tj^\o) 
and r^'*(0) with an insertion of the operator Ui and Tj^, respectively, i.e. 

where B and A are appropriate constant tensors such that Zu = = 1 at tree level. Then, 
we compute the RG functions rju and tjt defined by 



, . din Zu,T 



dlnm 



= /?„^^ + A^^, (A6) 

Uo,Vo 



where f3u and f3y are the /9-functions. 

We computed the functions T\j'j.{0) to six loops. The resulting six-loop series of rju^Tiu, v) 

are 

/ N 2m 1 2 12 + 2iV 4 2 2 V- C/ i r^rr^ 

^) = - sTiv - r + " WTW + WT^"" ^ ^ ? ' ^ ^ 

2m 2 12 + 2iV 4 
' ^ 8 + AT 3(8 + iV)2 9(iV + 8) ^ 

where and are normalized so that 

8 + M 3 , , , 

= — - — ZuUo, mv = -— Z^vq, Zu,v = 1 + v), (A9) 

457r iDTT 

and the coefficients efj an are reported in Tables X and XI respectively. Note that u 
and V correspond to u and v in Ref. [27]. The RG dimensions yu and yx are obtained by 
yu,T = '^ + flu,T — V-, where rju^x is the value obtained by resumming the corresponding series, 
evaluating it at w = m*, w = w*, where {u* ,v*) is the stable FP. 

In the case of the REIM, i.e., in the hmit A" — > 0, one has yu — He — 1/^^- Indeed, 
{cpfipk'^iy^^ and {Y.i'^'i'^k^iY^^ are both finite and nonvanishing for A" — > 0. Therefore, we 
have for A" ^ 

{NUi^pu^iY''' = -(E v]Vk^^^''' + 0{N) = -{NE<p,<piY''' + 0(Ar). (AlO) 
j 

Using the Monte Carlo estimate v — 0.683(3) (Ref. [19]), we obtain yu — 1.464(6). The series 
for TjT was already reported in Ref. [44]. Its analysis provided the estimate yx — 2.08(3). 

For A^ = 2 the stable FP of the cubic theory is the 0(2) FP, so y^ = yr = 1.766(3) 
[26]. For A^ > 3 the 0(A") FP is unstable and the RG trajectories fiow toward another 
FP characterized by a discrete cubic symmetry. The analysis of the series, using the same 
procedure reported in Ref. [27] , gives the estimates 

yu{N = 3) = 1.774(7), yriN = 3) = 1.800(2), (All) 
yu{N = 4) = 1.696(8), yriN = 4) = 1.874(3). (A12) 
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APPENDIX B: SOME RENORMALIZATION- GROUP IDENTITIES 



In this Appendix we prove relations (43). Morever, we show that, in the hmit N 
the RG functions do not depend on M ior v — 0. 

Let us consider the Hamiltonian for Vq — and rewrite 



0, 



exp 
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"4! 
exp 



ij,ab 

1 / , o ^ — \ o ^ — T Ox 1 



dXdpadua 



(Bl) 



where a (resp. i) runs from 1 to M (resp. N) and A, Po, and Uai are auxihary fields. Then, 
let us consider the n-point irreducible correlation function. We will show that it has the 
form 



(00 



a\%\...anln 



(B2) 



where Qa are group tensors (products of Kroncckcr deltas), the scalar factors do not 
depend on M, and the sum runs over all possible independent group tensors. In terms of 
the auxiliary fields, Feynman diagrams contain loops of 0-fields and n/2 open lines of 0- 
fields connected by the auxiliary-field lines. The group factor associated with each diagram 
is computed as follows. One assigns indices ai to and a propagators and indices a to p 
propagators, considers the product of the factors V (reported below) associated with each 
vertex, and sums over all assigned indices. In order to prove the M-independence we will 
show that, because of the Kroneckcr 5's appearing in the diagrams, none of these sums is 
effectively performed in the limit iV — 0, so that no factor of M can appear. The factors V 
are given by: 



V{\ (f)ai, (f)bj) 



^ab^ij I 



V{pcAaiAb3) = -^5abAj, 
y{(7ck-,4>ai-,4>bj) — -^^^^abc^ijk-, 



(B3) 



where ahc (resp. ijk) run from 1 to M (resp. N). First, note that all loops must contain 
at least a (T00 vertex, otherwise by summing over the indices of the fields appearing in 
the loop one obtains a factor of N . As a consequence, all loops give rise to a very simple 
effective vertex for the auxiliary fields: 



(A . . . Apai ■ ■ ■ Pa^CTfeiii ■ ■ ■ CTbmim) ~ <^ai...a„6i...6,„(^n...j„ 



(B4) 



where we have only written the dependence on the group indices. Then, given a diagram, 
let us consider the reduced diagram in which all 4> loops are replaced by the corresponding 
effective vertices. The (Tai propagators form several connected paths. It is easy to convince 
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oneself that each of these paths must end at an open hne, otherwise, by summing over the 
indices i associated with the a hnes, one obtains factors of N. As a consequence, all effective 
vertices are connected by a propagators to the open lines. Therefore, by summing over 
the indices associated with the (j) and a propagators one obtains expressions in which all 
remaining indices (those related to p propagators) are equal to external ones and are not 
summed over. We have thus proved that correlation functions expressed in terms of the 
bare parameters do not depend on M. Since Rmn is M independent for N 0, this result 
extends trivially to the RG functions expressed in terms of the renormalized couplings. 

To prove identities (43) we now exploit the M independence. For M — 1 the theory 
corresponds to the REIM model with couplings Uq + Wq and yo- Since Rmn — Rn for 
N ^ 0, {u + v)/m is a function of uq + wq- The result follows immediately. 
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TABLE II. The five-loop series of = mdu/dm for M = 3. Setting f5u = J2i=ol^u,h where I 
is the number of loops, we report I3u,l up to Z = 5. 



I Am 






—u 


1 


u'' + j^uv + ^uw + T^i;ui + ^uy 


2 




3 


0.389923 u* + 1.16913 + 0.621933 v'-' + 0.10714 uu^ + 0.64302 w + 1.1669 + 0.41177 ui;^ w + 0.02168 1;-^ «) 
+0.35637 + 0.394263 uvv? + 0.048431 + 0.0730773 uw^ + 0.0254122 vw^ + 0.857364 2/ + 0.912169 2/ 
+0.235717 uv'^y + 0.95032 v?wy + 0.597082 uvwy + 0.023813 v'^wy + 0.292309 uw'^y + 0.0327432 vw'^ y + 0.46739 j/^ 
+0.230387 u V y^ + 0.299507 u w + 0.00418781 vwy'^ + 0.090449 u y^ 


4 


-0.447316 m" - 1.83254 u*t) - 1.7466 - 0.575817 - 0.076344 md" - 1.0079 - 2 .83639 u-^ v w - 1.9044 «^ 
-0.380099 MD^ MI - 0.0145582 1)4 «) - 0.885278 u)^ _ i,74084.u^ v w'^ - 0.695073 wf^ ui^ - 0.038822 t)-^ «)2 _ 0.37419 ^,,3 
-0.451159 Mt)«i» - 0.0547698 - 0.06557 m to* - 0.0232153 - 1.34386 u'^ y - 2.56167 u^vy- 1.2668 u'-' v'^ y 
— u.zzoyio (1 u y z.ouui4: u ^ y — o.vju^i u, u Uu y — u.ofUiU u, u tfy — "J - yj^o loo u y — _l . 'lyu / u u w y — ±.11001 ll o w y 
-0.0520405 ?/ - 0.349706 Mw^y - 0.0459449 w)^ y - 1.2213 u^y^ - 1.23262 u'^ vy'^ - 0.32717 uv^y^ - 1.5635 «^ m) 
-0.888055 uvwy^ - 0.0135895 v^wy'^ - 0.550107 uup- y^ - 0.0186856 vv? y^ ~ 0.476234 v? y^ - 0.252712 uvy^ 
-0.335525 uwy^ - 0.00289791 vwy^ - 0.0754467 u y'^ 


5 


0.633855 u"^ + 3.32208 u'='!; + 4.6464 + 2.53461 + 0.702909 i;* + 0.0816358 ui)" + 1.82714 «; + 7.04207 t- w 
+7.40684 u^v'^w + 2.97466 w + 0.499013 uv* w + 0.0162882 w + 2.22077 w'^ + 6.4639 vw^ + 4.81333 
+1 . 15478 uv^w"^ + 0.0555485 1)* lu^ + ^ 45911 „3 ^3 ^ 3.06664 vw^ + 1.35727 uv^ + 0.0935349 + 0.522432 
+0.661835 uvw^ + 0.0871471 w'^ + 0.0804704 uw^ + 0.0286439 vw^ + 2.43619 y + 6.81471 2/ + 5.57615 y 
+2.06187 u'^v^y + 0.29933 uv'^y + 5.92205 u*wy + 12.266 u^vwy + 7.27261 u'^v'^wy + 1.54143 wy + 0.037268 v*wy 
+5.83645 u^w'^y + 8.62368 u'^vw'^y + 2.93576 wu^ w'^ y + 0.111143 w'^ y + 2.7863 y + 2.37933 uvw^ y 
+0.148969 v'^w^y + 0.53647 uw'^y + 0.0809503 vw'^ y + 3.14287 y^ + 5.2433 u^vy^ + 2.89085 y^ + 0.557283 u v'^ y^ 
+6.2083411''' m; 1/2 + 7.3735 m2 „ ^ + 2.24663 mv^ ^ + 0.034735 ui i/^ + 4.45664 ^,2 -i^ 2.95624 uv w'^ y^ 

+0.0805995 1)2 iu2 y2 + 1.14677 11 lo'' y^ + 0.0642921 dui^ i/^ + 1.93839 it"' + 2.12342 u^^y^ + 0.61093 ui;^ y^ 
+2.79532 «2 10 1/3 _^ 1.66738 11 f 101/'' + 0.0148483 id j/-' + 1.08348 uui^ yS + 0.020416 1; u)^ y3 ^ 0.630396^2 ^4 
+0.362335 uvy"^ + 0.489119 uwy* + 0.0022041 vwy^ + 0.0874933 u y^ 



TABLE III. The five-loop series of /3„ = mdv/dm for M = 3. Setting /3„ = J2i=oPv,h where I 
is the number of loops, we report up to Z = 5. 



3^ +1,2 + 2^ +-±nL 



77MDm Alv'^w 23 vw'-' 77 uvy 400 v'' y 23vu>y 92 vy'' 

2ifi 222 aafiX, ins flfl 2J,fi Saj ai Z22_ 

.,4 .1 1 01 «Q ....^ «„ _L 1 '7nn'7n 



0.916668 M''i, + 2.35093 w'' i)'' + 1.40888 uu^ + 0.28295 1;* + 1.2163 u'' i)to + 1.79079 uu'' «) + 0.542156 1;-" w + 0.49351 uuu)'' 
+0.367154 1,2 u,2 _|_ 0.0730773 vw^ + 1.62177 u^vy + 1.7731 u 1,2 + 0.508415 y + 1.31602 uvwy + 0.746208 v^wy 
+0.292309 vw'^y + 0.63527 it vy'^ + 0-346201 1,2 y2 _^ 0.299507 vwy'^ + 0.090449 v y^ 

- 1.22868 u*v- 4.3522 v'-' - 4.33525 - 1.75357 u - 0.270333 1," - 2.26103 tt^ i, lu - 5.40209 tt^ i,^ w - 3.5142 u i;3 « 
-0.712012 v'^w- 1.53326 vw'^ - 2.40612 uv^ - 0.798138 up' - 0.486101 uvuP - 0.401462 1;2 - 0.06557 v 
-3.0147 v?vy- 5.65057 1*2 „2 ^ _ 3 27089 ii i;^ y - 0.650693 v'^ y - 4.0887 u^vwy- 5.03042 11 1;2 «; 1/ - 1.47787 v^wy 
- 1.9444 lii; 11,2 y - 1.26296 1;2 w'^ y - 0.349706 1; ii,^ y - 2.1004 ii2 i,i/2 _ 2.31802 iii;2 y2 _ 0.671001 1;^ y2 _ 2.02679 uvwy^ 

■ 0.616942 11 1, y3 - 0.350699 1,2 - 0.335525 1, u, y^ _ 0.0754467 1, y* 



1.97599 11-'' 1, + 9.0062911*1,^ + 12.7239 if' 1,^ + 7.96583 «^ t," + 2.46634 ut--'' +0.312556 1," + 4.6941711" 1,11, + 15.6644 u'' 1,2 11, 
+16.0384 ii2 i;3ii, + 6.69086 11 1--* lu + 1.04399 1,'^ 11, + 4.50574 1, ui2 + 11.0238 ii2 1,2 «;2 + 7.53313 « y-' iy2 + 1.60193 1,* iu2 
+2.28761 ii2i;ii,=' + 3.82594 « ^2 iu-' + 1.35045 1,^ ii,^ + 0.642513 11 1, 10* + 0.562428 1,2 ii," + 0.0804704 1, «,■"'' +6.25889 u'^vy 
+16.8745 11^ i;2 y + 15.2053 ii2 v'^ y + 6. 1052 u y + 0.947895 1,^ y + 12.0153 ii^ 1, 10 y + 23.5313 ii2 1,2 «, y + 14.2714 u i,^ u, y 
+2.9145 1,-* i«y + 9.15043 u2 1,11,2 y ^ 12.2315 11 1,2 11,2 y + 3.84191,-'' 11,2 y + 3.42673 111; iw^ y + 2.40324 1,2 ii,3 y + 0.53647 1,11;'' y 
+6.36312 ii3i;y2 + 11.045 ii2 1,2 y2 +6.4832 11 1,^ y2 + 1.31809 1,* y2 + 9.71176 ii2 „^j/2 _,_ 11.5283 u i;2 luy 2 +3.48488 i,3ii,y2 
+5.47297 ui;«,2j^2 +3.4119^2 ^2 y2 _,_ 1.14577 ^«,3 y2 + 3.01986 1*2 i;y3 _|_ 3 47754^^2 yZ _,_ 1,04671 1;3 y3 + 3.42367 u i; «, y^ 

+2.05758 v^wy^ + 1.08348 d «,2 j/3 + 0.771673 w d y"* + 0.461776 1;2 y'^ + 0.4891 19 vwy'^ ^ 0.0874933 v y^ 
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TABLE IV. The five-loop series of = m,dw/dm for M = 3. Setting fiy^ = ^i=Ql3w,h where 
I is the number of loops, we report (5yj^i up to i = 5. 



:iuw I Avw I ,..2 I '^wy 

2 11 2_ 



223 uvm _ 244 w _ 131 u _ 47 v w' _ 95 _ 77u-wy _ 184 v w y _ 50 w'' y _ 92 ni y'' 



21£ 222 22S2 IQS 21fi , a: ., aj Z22 „ ^ 

0.916668 w + 1.37496 '10 + 0.362132 md^ w + 0.042103 1)^ w + 1.98317 + 1.74277 mvw^ + 0.271146 

+1.48661 ifU)-^ + 0.676067 uw)^ + 0.389923 w'^ + l.mm wy + 0.939981 m d w) j/ + 0.164277^)2 + 2.27978 mw^j/ 
+0.732074 y _|_ 0.857364 ?/ + 0.63527 uwjy^ + 0.217823 vwj/^ + 0.467389 u)^ j/^ + 0.090449 wj/^ 



-1.22868 it* «) - 2.66348 -u'^t) a- - 1.4268 «^ v'' w - 0.347251 m •u'^ ty - 0.0326696 v'^ w - 3.58788 - 5.25808 t; 

-1.93572 m (L-^ - 0.244281 'f'"' u.^ - 4.09897 li)-^ - 3.9245 u f u)'' - 0.759482 ?)2 - 2.17101 ww'' - 1.07631 
-0.447316 «)5 _ 3.0147u=^ u) 2/ -3.33336 u^uwj y - 1. 18975 u i)^ u,^ _ 0.154387 1)^ _ 6.5664 ti,2 ^ _ 4.83085 ui)«)2 

-0.899584 v'^w'^y- 5.02574 uw^y- 1.95255 d ii)'' y - 1.34386 lu* 2/ 

-2.1004 u2tj,y2 _ 1.46007m I) _ 0.2864021)2^4,^2 _ 3.1138 uju^ _ 1.15025 1)«)2 j,2 _ 1.2213 u)3 y2 _ 0.616942 wwy^ 

-0.244019 vwy^ - 0.476234 w'^ y^ - 0.0754467 w 1/^ 

1.97599 u''w + 5.70624 d tu + 4.72927 w + 1.77266 w + 0.357776 uv*w + 0.0327712 i)" w + 7.28706 u* w'^ 
+15.3337^3 I) u;2 + 9.42002^2 ^2 ^2 ^ 2.47241 ud^ uj^ + 0.26882 1)^ u)^ ^ 11.2179^2 + 17.1513 ■yju^ + 7.14394 ui)^ 
+0.985494 1)3 «)3 + 8.985399 u'^ + 9.21597 uvw'^ + 1.99704 1;^ ^4 _,_ 3 72266 uw^ + 1.9646 d lo^ + 0.633855 
+6.25889 u'^wy + 10.5195 u^vwy + 5.91865 ^2 ^2 ^ ^ _,_ ;^ g;^48 uv'^wy + 0.18753 v'^wy + 18.525 w'^ y + 22.759 v? vw'^ y 
+8.9157214 1)2 w^y + 1.29356 1)=* w'^ y + 21.5756 1*2 nfl y + 18.0042 « d lo'^ y + 3.70671)2 ^3 _,_ 11.6445 iiu;"' 1; + 5.06192 uiu^y 
+2.4361911)5 1/ + 6.3631211-^ lu 1/2 4. 7.15909 u^vwy^ + 2.9287 11 1>2 lu i/2 + 0.453078 1)^101/2 + 14.4006 ii2 id2 y'^ 
+11.226111)11)2 1/2 + 2.42873 1)2 ii;2 i/2 ^ 11.4247 u ii)^ y^ + 4.69358 1) ii;^ i;2 + 3.14287 id"' i/2 + 3.01986 u^wy'' 
+0.5663871)2 ID + 4.73169 ««)2y3 + 2.0427 dio2 j^3 _,_ 1.93839 lo^ j^s _,_ 0.771573 uio i/* + 0.355723d 10 + 0.630396 iu2 j/" 
+2.45617 uvwy^ + 0.0874933 id y^ 



TABLE V. The five-loop scries of (3y = mdy/dm for M = 3. Setting Py = J2i=q (^y,h where I is 
the number of loops, we report fiy^i up to / = 5. 



9tim I Suy , Vivy , A w y , 2 
_LJ Z 2 Z U Z 2_ "^^ w 



9uvw 6'.iv^w 27 V w'^ 18.5 it'' y 439 u ^' y 17.561?'')/ 185 u ^) j/ 185^11'^ 

11 „ 121 „ 44 „ 216 „ 297 3267 108 99 

185 U) y _ WAuy-' _ 832 v y^ _ 104 iii y^ _ 308_yj^ 
891 m „ 729 



1.25483 ii2 DID + 1.62749 •ud2 id + 0.44594 d'^ id + 1.84622 it d id2 + 1.1317^2 + 0.690211 d id'' + 0.916668 ii''' y 
+2.49036 ii2 Dy + 1.80879 11 d2 i/ + 0.438494 v'^ y + 2.75 ii2iDy + 5.96358 «d id y + 2.36751 d2 wy + 2.75 uw'^ y + 3.2512 vw^y 
+0.91667iD3y + 2.133 ii2 + 3.1195 uvy'^ + 1.13436 d2 y2 + 4.26599 11 id i/2 + 3..38291 d id y2 + 2.133 w'^ y^ + 1.47806 ui/^ 
+1.07495 py^ + 1.47806 id i/-^ + 0.351071/" 



2.17122 DID -4.35781 u2 d^ id - 2.48256 it d^ id - 0.48565 d* id - 4.77724 1(2 did2 - 6.00211 ud2 id2 - 1.72925 d^ id^ 
-3.62343 UDiD^ - 2.2762 d2 11)3 - 0.951639 did" - 1.22868 w" y - 4.59345 dj/ - 5.3004 d^ y - 2.55397 MD^y 
-0.464358 d" y - 4.91474 ii^ id y - 16.29971*2 did y - 13.3566 11 d2 i« y - 3.38345 d^ wy- 7. 37211 v? y - 17.4996 11 d ii)2 y 
-7.45225 d2 11)2 y _ 4.91474 u u)-^ y _ 6.09451 d id^ y - 1.22868 id" y - 3.89927 Ii3y2 _ 9.02889 li^ Dy2 - 6.52141 iid2 y2 
-1.58095 D3y2 - 11.6978 u2 id y2 - 19.0578 UD id y^ - 7.19018 d2 id y2 - 11.6978 wid^ y2 - 9.94972 d id2 y2 - 3.89927 id^ y2 
-4.23696u2y3 -6.1212iiDy3 -2.22589d2 y3 - 8.47392 w id y^ - 6.37921 DiDy^ -4.2369611)2 yS 

-2.03309My" - 1.47861 py" - 2.03309 Wj/" - 0.376527 y^ ^ 

4.2429179 u" d id + 11.7357 d^ id + 10.5256 11^ d^ id + 4.10472 11 d* id + 0.618728 d" id + 12.4807 d id^ + 24.3706 11^ v'^ w'^ 
+14.3671 11 d3 11)2 + 2.82095 v* 'ufl + 14.2955 ii2 d lo^ + 18.2581 11 d2 lu^ + 5.42606 d^ id^ + 7.5579084 11 d id" + 4.86798 id" 
+1.55088 did"' + 1.97599 ii^ y + 9.47772 u" d y + 15.161 ii'^ d2 y + 11.1288 d^ y + 4.00642 11 d" y + 0.582751 d"' y 
+9.87996 11" ID y + 44.5209 Dioy + 56.6807 ii2 d2 iDy + 28.7278 ud^ id y + 5.,3291d" wy ^ 19.760 ii^ id2 y + 71.398 u^vw^y 
+61.8264 11 d2 11)2 y + 16.0162 d^ id2 y + 19.760 ii2 id^ y + 49.2653 md id^ y + 21.8335 d2 id'^ y + 9.87996 mid" y + 12.6315 d id" y 
+1.97599iD5y + 7.98749 11" y2 + 25.6071 „ y2 + 28.2134 ii2 „2 y2 _,_ 13.5424 ud^ y2 + 2.46225 d" y2 + 31.9499 ^Dy^ 
+80.4247 ii2 D 10 y2 + 60.7793 uv^wy^ + 14.9778 d^ 10 y2 + 47.9249 ii2 ^2 y2 _,_ 82.7197 11 d id2 y2 + 32.0629 d2 id2 y2 
+31.9499 11 id3 y2 + 28.2945 d 11)3 y2 + 7.98749 id" y2 + 11.9097 ii3 y3 + 26.2915 ii2 d y3 + 18.9292 iid2 y^ + 4.58889 d3 y3 
+35.729 1(2 ^y3 + 53.9545 MD ID y3 + 19.9234 d2 id y3 + 35.729 wid2 y3 + 27.6713 did2 y3 + 11.9097 id3 y3 +8.72598 1*2 y4 
+12.5648 11 Dy" + 4.56903 d2 y" + 17.452 mid y" + 12.881 did y" + 8.72598 id2 y" + 3.24652 wy^ + 2.3611 Dy^ + 3.24652 id y^ 
+0.495548 y^ 
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TABLE VI. The five-loop series of Pu = mdu/dm for M = 2. Setting /?„ = Y.i=o Pu.i, where I 
is uumlicr of loops, we report ;^„j up lo / = 5. 



95 20ii^i> 92 uv''' 'Zbu'-'m Aluvw 2 v'^ w 23 u w'^ v w'^ 50 u'' y 92 u v y 23. uw y 92 u y'' 



7B 21 31) aj i!JS aj Z22_ 



0.389923 m"' + 1.02884 d + 0.553406 + 0.101002 + 0.643023 «^ «) + 1.22281 -uw + 0.444702 mv^ «;+ 
0.025201 v-^ w + 0.35637 + 0.429803 m d lo^ + 0.0577681 v'^ w'^ + 0.0730773 u w'-^ + 0.0279534 vw'-^ + 0.857364 y 

+0.922343 V y + 0.252504 uv^ y + 0.95032 u^wy + 0.646426 uvwy + 0.028814 v^wy + 0.292309 uw'^ y + 0.03602 vw'^ y 
+ 0.248819 «Dj/^ + 0.299507 u '11) + 0.00460659 u 'U) + 0.090449 « j;'^ 



-0.447316 m" - 1.61263 m"'?) - 1.50485 m'^ t)^ - 0.539552 v'-^ - 0.0779026 Mt)" - 1.0079 - 2.88772 ?t-^ w; 
-2.0261 M^ 1)2 ui - 0.444498 Mf'^ m) - 0.0187035 i)"* m) - 0.885278 m* m)^ - 1.88177 m^ i)m)2 - 0.82283 md^ u)^ - 0.051806 1)^ u)^ 
-0. , 374191 M^ M)3 - 0.493586 M t) M)3 - 0.0658781 m^ m)'^ - 0.06557 mm)"' - 0.0255368 m m)"* - 1.34386 m"' ?/ - 2.50809 m^ My 
-1.34888 M^M^ y - 0.259675 mm* y - 2.36074 m* m)j/ - 3.28057 m^ mm; j/ - 1.0336 mm^ m; i/ - 0.031355 m* m) y - 1.49676 m^ m)^ j/ 
-1.2471 m«;u)2 2/ -0.06253 m2m;2,; - 0.349706 mm)3 _ 0.0505394 m m)^ _ i.2213u3 j/2 _ l.32086u2My2 _ 0.381274 mu^ ^^2 
-1.5635 u^tu 2^2 _ 0.966102 u MM) 1/2 _ 0.0164433 m^ wj/^ - 0.550107 m m)^ _ 0.0205542 m m)^ _ 0.476234 y3 
-0.274796 MMy^ - 0.335525 mm) y* - 0.0031877 m w - 0.0754467 My 



0.633855 u** + 2.92343 m'^m + 3.946 u-*M^ + 2.28569 m^ + 0.683062 m"* + 0.0863513 mm'' + 1.82714 u'' m; + 7.05664m* mm) 
+7.67985 u^v'^w-\- 3.31184 w -\- 0.600884 m m* m) + 0.021403 w + 2.22077 m* m)^ + 6.8986 m^ m to^ + 5.54476 v? 
+1.44276 u)2 _,_ 0.0762973 m* + 1.45911 + 3.34723 vw^ + 1.6199 um^ + 0.122458 + 0.522432 

+0.725783 M M m;* + 0.104985 m)* + 0.0804704 u m;^ + 0.0315083 vw^ + 2.43619 y + 6.57667 m y + 5.71423 y 
+2.27687 u'^v^y + 0.359797 m m* y + 5.92205 m* m) j/ + 12.9278 m* m m; y + 8.30663 v?v'^wy + 1.91539 uv^wy 
+0.0515184 v*wy + 5.83645 y + 9.38178 m^ m m)2 y + 3.48908 uv^w^y + 0.146069 w'^ y + 2.78631 m^ y 

+2.60534 M M M)3 y + 0. 17948 m^ m)* y + 0.53647 m m)* y + 0.089045 m m)* y + 3. 14287 m* y^ + 5.50064 m* m y^ + 3.27895 m^ m^ y^ 
+0.689359 M M* y2 + 6.20834 m* m) y^ + 8.00363 m^ m m) y^ + 2.6647 m m^ «> y2 + 0.0459087 m* m; y^ + 4.45664 m^ y^ 
+3.23343 MM M)2y2 + 0.0973883 m2 m)^ y2 + 1.14677 uw'-^ y^ + 0.0707213 mm)* y^ + 1.93839 m* y* + 2.30259 m^ My* 
+0.723961 mm2 y* + 2.79532 m^ M)y* + 1.82321 mm M)y* + 0.0179664 m^ m) y* + 1.08348 m m)2 y* + 0.022458 M u)^ y* 
+0.630396 m2 y4 + 0.396144 mm y"' + 0.489119 m m) y* + 0.00242451 vwy'^ + 0.0874933 My^ 



TABLE VIL The five-loop series of = mdv/dm for M = 2. Setting = Ez=o /^M,i, where Z 
is the number of loops, we report Py^i up to Z = 5. 



/ 


!3v.l 





— V 


1 


— 5 T 1 2 

M + ^MM + -VW + -My 


2 


185 u^u 181 UD^ 272 77 u w to 47d^to 23uto^ 77 uvy 40v''y 23 v w y 92 vy'' 
21 fi 1.15 fi7.^ ms Sn 21 fi SI Rl Rl 723 


3 


0.916668 M*M + 2.25286 m'' m^ + 1.44833 mm* + 0.314917m'' + 1.21633 u'' mm) + 1.9382 mm'' m) + 0.635652 m* m) + 0.49351 mm m)'' 
+0.401926 m2 m)2 + 0.0730773 m m)* + 1.62177 m^ m y + 1.90822 uv^ y + 0.593025 m* y + 1.31602 mm m) y + 0.815647 v'^wy 
+0.292309 M M)2 y + 0.63527 uvy^ + 0.378517 m^ y^ + 0.299507 M M) y^ + 0.090449 m y* 


4 


-1.22868 M*M -4.14839 u*M^ - 4.30604 u^m* - 1.88967 u m* - 0.317928 m'' - 2.26103 m* mm; - 5.70167 m^ m^ m) - 4.000 mm* m) 
-0.885211 M* M) - 1.53326 m^ m m)^ - 2.62769 m m^ m)^ - 0.952554 m* m)^ - 0.486101 m m m;* - 0.440264 m^ m)* - 0.06557 v 

-3.0147 vy - 5.8948 m^ m^ y - 3.71358 mm* y - 0.8076.36 m"* y - 4.0887 m^ mm) y - 5.48266 mm^ M)y - 1.7649 m* M)y 
-1.9444 MM M)2 y - 1. , 38387 m2 w'^ y - 0.349706 m m)* y - 2. 1004 m^ My^ - 2.53023 mm^ y^ - 0.801818 m* y^ - 2.02679 m m m) y^ 
-1.27603 m2 M)y2 - 0.550107 MM'2y2 _ 0.616942 m My* - 0.384175 m^ y* - 0.335525 mm; y* - 0.0754467 m y"* 


5 


1.97599 m='m + 8.54783 M^M^ + 12. ,3856 m* m* + 8.36256 m^ m" + 2.82533 mm" + 0.391102 m'' + 4.69417 m"' mm) + 16.29.55 m* m^ m) 
+17.8673 m2 m*m) + 8.13124mm4m; + 1.38643 m-''' m) + 4.50574 m* m m)2 + 11.9507 m^ m^ m)2 +8.92856 mm* w'^ + 2.07519 m* m)2 
+2.28761 m2 MM)* + 4.19835 mm^m)* + 1.62409 m*m)* + 0.642513 mm m)"' + 0.617553 m^ m)"' + 0.0804704 m m)" +6.25889 M"'My 
+17.3148 M* m2 y + 16.873 m^ m* y + 7.4013 u m"' y + 1.25617 m^ y + 12.0153 m* m m; y + 25.4165 m^ m^ m) y + 16.8835 m m* m) y 
+3.77286 m"' M)y + 9.15043 M2MM)2y + 13.414 mm^ m)^ y + 4.61743 m* m;^ y + 3.42673 m mm;* y + 2.6376 mi* y + 0.53647 m m)"' y 
+6.36312 M*My2 + 11.940 M^ u^y^ + 7.66554 mm* y^ + 1.70545 m'' y^ + 9.71176 m^ m m) y^ + 12.6,39 m m^ -u; y^ +4.18792 V^wy^ 
+5.47297 MMM)2y2 + 3.74384 m2 y^ + 1.14677 m m;* y^ + 3.01986 My* + 3.81205 mm^ y* + 1.2579m* y* + 3.42367 u m m) y* 
+2.25788 m2 m) y* + 1.08348 m m)^ y* + 0.771673 m m y"" + 0.506741 m^ y* + 0.489119 m m) y* + 0.0874933 m y^ 
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TABLE VIII. The five-loop series of = mdw/dm for M = 2. Setting = Y^i^q Pw,l, where 
I is the number of loops, we report ^^,1 up to Z = 5. 



77 u w y 92vwy 50 y 92 w y 

2ifi „ 27 fijs ins .fln -^rnr _ si 4nR aj I22_ 



0.916668 u-^w + 1.33266 u'' vw + 0.406785 ui)^ w + 0.0505995 w + 1.98317 w'^ + 1.90301 uvw'^ + 0.324868 v'^ w'^ 
+1.48661 utu^ + 0.743674 t)tt)^ + 0.389923 + 1.62177 wy + 1.04026 u u «; j/ + 0.194876 1;^ w 2/ + 2.27978 wtt)^ y 
+0.805282 vw'^y + 0.857364 y + 0.63527 uwy^ + 0.239606 vwy"^ + 0.467389 y'^ + 0.090449 w 



-1.22868 u'^w- 2.55616 u-^ vw - 1.47403 v'' w - 0.384865 uv' w - 0.0404956 v'' w - 3.58788 - 5.59606 v w'-' 
-2.25333 11)2 - 0.310265 - 4.09897 ,^,3 _ 4,s3198uvw^ - 0.921371 v'^ - 2.17101 um)" - 1.18394 iJio". 

0.447316 lu'^ - 3.0147 wy - 3.53479 u'^ vw y - 1.37253 m ui j/ - 0.196966 1)^ - 6.5664 ,„2 _ 5 .33631 u v w'^ y 
-1.09183?)2 w'^y- 5.02574 MW)'^ 2/ - 2.1478««)'^ y - 1.34386 tu"*?/ - 2.1004 «)j/2 _ i,Qii2uvwy'^ - 0.348387 v'^ w y^ 
-3.1138 ?tw)2j^2 _ 1.26528 ?)W)2j^2 _ i,2213 w^y'^ - 0.616942 m «i j/^ - 0.26842 tjwj/^ - 0.476234 «)2 _ 0.0754467 y"* 



1.97599 u^w + 5.43635 u'^vw + 4.70998 u'-^ v^ w + 1.94106 w + 0.434203 uv'^w + 0.0435452 v^ w + 7.28706 w'^ 
+16.0529 M^t)W)2 + 10.7128 «2 „2 ^2 + 3.09443 t«2 _,_ 0.3684 ui^ _,_ 11.2180^3 u)^ + 18.7355 _,_ g g]^4]^;^ „^2 ^3 

+1.30212 ?)3 w^ + 8.9854 m2,„4 + 10.1884mi;w4 + 2.4288 1)^ ty4 + 3.72266 mw)^ + 2.16106 + 0.63385 + 6.25889 
+ 10.9184 u^vwy + 6.74643 u^v^wy + 2.0261 uv^wy + 0.256761)" wy + 18.5254 m3 ^2 _|_ 24.8502 u^vw^y 
+10.7457 ui;2 w'^ y + 1.7078 v"^ w'^ y + 21.5756 w"^ y + 19.9199 u v w'-'' y + 4.50834 w'^^ y + 11.6445 uw" y + 5.56811 vw'^y 
+2.4362 6.3631 u3iuy2 _|_ 7.8398 „t„ ^2 +3.53015^1,2 w y'^ + 0.59754 v^ w ,y^ + 14.401 u'^ w^ y'^ + 12.417 u 1; «;2 ^^2 

+2.95262 1;2 ^2 ^2 + 11.4247 uu)^ 2/2 _,_ 5 ^6294 1; u)^ + 3.14287 u;" 2/2 _,_ 3 01986 u2 „,j,3 _,_ 2.7i4i7u„^ j^3 ^ Q g8gQ3 ^2 ^ ^3 
+4.73169 u 14)2 y3 + 2.24697 D«)2j^3 ^ 1.93839 w3j/3 + 0.771673 uiuj/^ + 0.391295 i^iuy* + 0.630396 «)2 ^ q 0374933 ^^^5 



TABLE IX. The five-loop series of Py = mdy/dm for M = 2. Setting /3„ = J2i=oPy,l^ where I 
is the number of loops, we report Py^i up to Z = 5. 



-y 

in 2 Z 5 1 2 ~^ » 

9 uv u} 3 V w 27 V w 185 u y 208 uvy 416 V y 185 uwy 37 v w y 

10 ^ 5 o 40 .,216 ^ 135 , 675 108 18 

185 y _ 104 Mi;-^ _ 416 _ 104 mi;-' _ 308 y'' 

21 .) aj 405 a: _J22 „ „ — „ , , 

1.38031 u''vw + 1.8544 wu^ w + 0.559618 w + 2.03085 uvw'' + 1.36855 v'' + 0.759232 vw'' + 0.916668 w3 y 

+2.5596 v?vy + 2.05514 uv'^ y + 0.548037 v^ y + 2.75u'^ wy + 6.52715 uvwy + 2.86347 ?;2 w j/ + 2.75 u 1^2 y + 3.57632 vw'^y 

+0.916668 2/ + 2.133 u2 y'^ + 3.42523 uvy'^ +1.370091)2 j/2 + 4.26599 uiu 2/2 + 3.7212 1; lu 2/2 

+2.133 11)2 y2 ^ 1.47806 M2/^ + 1.18245 1;2/^ + 1.47806^2/^ + 0.351072/^ 



-2.38834 m3„«, -4.91662 i;^ lu - 3.05722 mi;3 u, - 0.652482 1;* lo - 5.25496 i; id'-" - 7. 19501 m d'-" ii)'-" - 2.27241)3 m,^ 
-3.98577 111; 10^ - 2.75879 1)2 11,=* - 1.04681)111" - 1.22868 11" 1/ - 4.67913 d 1/ - 5.84436 ii2 1)2 y - 3.1024 11 1''^ 1/ 
-4.91474 11^ ID 2/ - 17.653 ii2 1)11,1/ - 15.9343111)2 11,1/ - 4.43201 id 1/ - 7.37211 ii2 id2 1/ - 19.2778 iid w)2 y _ 9.03619 d2 id2 y 
-4.91474 11 iD^ 1/ - 6.70396 vw^y- 1.22868 id* 1/ - 3.89927 i/2 - 9.7406 u^vy^- 7.756 md2 i/2 - 2.06827 d^ i/2 
-11.6978 ii2 ID 1/2 - 21.0095 11 DID 1/2 - 8.71866 ^2 id 2/^ - 11.6978 uid2 2/2 - 10.9447 vw^ y'^ - 3.89927 id^ i/2 - 4.23696 «2 y^ 
-6.74859 11 D 1/3 - 2.69944^2 y"' - 8.47392 u lu i/^ - 7.01714 1; id 2/* - 4.23696 id2 yS _ 2.033091*2/"* - 1.62647 1; 2/" 

-2.03309 ID 1/" - 0.376527 1/"' - 0.62048 d"i/ ^ 

4.6672 11" D w) + 13.1278 ii3 D^ lu + 12.6693 d3 id + 5.41225 iid" id + 0.895128 d» id + 13.7288 ii'^did^ + 28.8549 11^ d^ id^ 
+18.6536 11 D'^ iy2 +4.02044 d" 11,2 + 15.7251 ii2 d lo^ + 22.1012 u d2 id'^ + 7.22562 d'^ ly'^ + 8.3137 11 d id" + 5.90266 d2 id" 
+1.70597 did"' + 1.97599 + 9.58498 11" d 1/ + 16.3791 d2 1/ + 13.2026 ii2 d'^ 1/ + 5.24509 iid" y + 0.839214 d^ 1/ 
+9.87996 u" ID 2/ + 47.8346 u^vwy + 66.4888 ii2 d2 id y + 37. 149 u d^ id 2/ + 7.5788 d" w y + 19. 7599 «,2 + 78.284 ii2 d id2 y 
+74.7963 ii 1)2 u)2 2^ + 21.317 1^2 y + 19.7599 ii2 lu^ y + 54.3088 uvw^y + 26.4953 d2 u,3 y _,_ 9.37996 u id" 1/ + 13.8947 d id" 2/ 
+1.97599 lu^ 2/ + 7.98749 ii" 2/^ + 27.2985 li^i; 2/2 + 33.00661*2 ^2 ^2 ^ 17.4336 ui;3 y2 + 3.49673 d" 2/^ + 31.9499 u, 2/2 
+88.09741*2 did 2/^ + 73.5089 ii?;2 t„ 2/^ + 19.9292 id 2/^ + 47.9249 1*2 k,2 2/2 + 91.2671 iiDiw2 2/^ + 38.9216 1;2 iu2 j/2 
+31.9499 liiD^ 2/2 + 31.124 DiD^ 2/^ + 7.98749 lu" 2/^ + 11.9097 li^ 2/^ + 28.80821*2^2/^+22.89051*1)2 j/3 ^0414^3 2/^ 
+35.729 1*2 iu 2/^ + 59.5478 lii; ID + 24.1885 v'^wy^ + 35.729 i*id2 y^ + 30.4384-!; id2 y^ + 11.9097 lu^ 2/^ + 8.72598 1*2 2/" 
+13.8681 1*1; 2/" + 5.54723 1;2 2/" + 17.452 1* id 2/" + 14.1691 did 2/" + 8.72598 id2 2/" + 3.24652 1(2/5 + 2.59721 1; 2/^ + 3.24652 id 2/^ 
+0.495548 2/^ 
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TABLE X. The coefficients eYj, cf. Eq. (A7). 




(AT + sfeY^ 


3,0 
2,1 
1,2 
0,3 


-18.3128 - 3.43328 + 0.216746 N'-' 

-9.15642 - 0.17027 Af 

-1.17334 

-0.0443103 


4,0 
3,1 
2,2 
1,3 
0,4 


140.799 + 37.5734 N + 1.03627 N'^ + 0.0943426 

93.8662 + 5.52597 N - 0.0781363 A^^ 

18.4511 - 0.0667897 AT 

1.40677 

0.0395196 


5,0 
4,1 
3,2 
2,3 
1,4 
0,5 


-1340.07 -416.717 AT - 17.6226 Af^ + 0.911281 Af^ + 0.0508337 AT" 

-1116.73 - 98.6844 AT + 1.52723 Af^ - 0.0301952 

-298.289 - 2.54766 AT - 0.0492195 AT^ 

-35.2065 + 0.140508 AT 

-1.98589 

-0.0444004 


6,0 
5,1 
4,2 
3,3 
2,4 
1,5 
0,6 


15651.3 + 5665.65 N + 433.687 A^-^ + 1.06755 A^'^ + 0.679106 Af" + 0.031393 AT'' 

15651.3 + 1935.63 N + 8.74297 N'^ + 0.581411 N'-^ - 0.00927903 W 

5294.38 + 134.776 N - 1.13059 A^^ - 0.038664 AT^ 

848.418 - 2.51108 A^ + 0.0323227 N"^ 

72.4799 - 0.318348 AT 

3.24291 

0.0603632 



TABLE XL The coefficients e^-, cf. Eq. (A8). 






2,0 

1,1 
0,2 


-18.3128 - 3.43328 AT + 0.216746 AT^ 
-3.09273 + 0.216746 AT 
-0.0337239 


3,0 
2,1 
1,2 
0,3 


140.799 + 37.5734 AT + 1.03627 AT^ + 0.0943426 AT-^ 
39.0459 + 1.53797 N + 0.12579 AT^ 
2.66843 + 0.0769829 N 
0.0716893 


4,0 
3,1 
2,2 
1,3 
0,4 


-1340.07 - 416.717 A^ - 17.6226 N'^ + 0.911281 Af'^' + 0.0508337 AT" 

-497.159 - 32.4255 AT + 1.57919 Af2 +0.0847229 

-53.6464 + 0.443225 N + 0.062623 N'^ 

-2.39176 + 0.0256721 N 

-0.0421612 


5,0 
1,1 
3,2 
2,3 
1,4 
0,5 


15651.3 + 5665.65 N + 433.687 A^^ + 1.06755 Af'^' + 0.679106 AT* + 0.031393 N'' 

7460.04 + 849.888 N + 0.972279 N'^ + 1.37677 + 0.062786 AT" 

1175.99 + 25.2714 AT + 1.12305 Af^ + 0.0567755 

88.2226 + 0.60426 N + 0.0297911 

3.53359 + 0.0136874 AT 

0.0607723 
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